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We study a solid plate plunging into or being withdrawn from a liquid bath, to 
highlight the fundamental difference between the local behavior of an advancing or 
a receding contact line, respectively. It is assumed that the liquid partially wets the 
solid, making a finite contact angle in equilibrium. In our hydrodynamic description 
which neglects the presence of the outer gas atmosphere, an advancing dynamic 
wetting line persists to arbitrarily high speeds. The receding wetting line, on the 
other hand, vanishes at a critical speed set by the competition between viscous and 
surface tension forces. In the advancing case, we apply existing matching techniques 
to the plunging plate geometry, to significantly improve on existing theories. For the 
receding contact line, we demonstrate for the first time how the local contact line 
solution can be matched to the far-field meniscus. In doing so, we confirm our very 
recent criterion for the instability of the receding contact line. The results of both 
the advancing and the receding cases are tested against simulations of the full model 
equations. 

I. INTRODUCTION 

A number of recent experiments have tested the stability of forced advancing 

and receding contact lines under conditions of partial wetting. For example, if a solid plate 
or fiber is plunging into a liquid bath to be coated (advancing contact line), the speed can 
be quite high (m/s) 4], while maintaining a stationary contact line. In the opposite case of 
withdrawal (receding contact line) j2l], a stationary contact line is observed only for very low 
speeds, and a macroscopic film is deposited typically at a speed of only a few cm/s. 

The description of a moving contact line is complicated by the fact that the Navier- 
Stokes equation with standard no-slip boundary conditions Q breaks down near it, because 
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FIG. 1: A schematic of the setup: a plate is being withdrawn at an angle 9 on the right, and 
being pushed into the fluid on the left. Only the flow to the left of the plate is considered, which 
is assumed long enough for the coupling to the other side to be ignored. At the contact line 
the microscopic slope of the interface is /i'(0) = Oe- Since the interface is highly curved, this 
microscopic contact angle is only observed on a scale of nanometers (expanded region on the left). 
The depression of the advancing contact line (left) relative to the level of the liquid bath is A. The 
contact line is stable at any speed, while for the receding contact line (right) there exists a critical 
speed Uc above which the contact line vanishes. Instead, the plate is covered by a thin film. 

such a hypothetical flow would produce an infinite energy dissipation 1?^. Instead, some 
microscopic length scale A must be invoked that cuts off this singularity, which in this paper 
we are going to take as a slip length. As a result, the local flow near the contact line is 
characterized by a typical length of about a nanometer [3], which has to be matched 
to the macroscopic flow away from the solid. In this paper, we apply a matching method 
developed for spreading drops [10] to the plunging plate, and test the result by comparing 
to numerical simulations. This very signiflcantly improves the results of earlier calculations 
for the same problem . For the opposite case of a receding contact angle we flnd that 

a new matching procedure is needed. The results conflrm our very recent criterion for 
the instability of a receding contact line. 

Figure n illustrates the geometry to be considered in the present paper. On the left, a 
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solid plate is pushed into a pool of viscous liquid. As a result, the interface deforms and 
the contact line is pushed downward relative to its equilibrium position. Within the present 
model, this advancing contact line is stable at any speed. If the plate is withdrawn from 
the bath, fluid is pulled up with the plate, and a new contact line position is established. 
However, this state is realizable only at speeds below a critical speed Uc- For U > Uc 
the contact line is no longer sustainable and continues to move up the plate Q|. In the 
stationary state, the plate is covered by a thin film, first described by Landau, Levich, and 
Derjaguin (LLD) HQ- 

To understand the fundamental difference between pushing a 
plate and pulling it out, one has to consider the matching between the region very close to 
the contact line, and the capillary profile away from it. Qualitatively, the difference in the 
behavior of the advancing and of the receding contact line makes sense: if the plate is pushed 
into the liquid, the interface is bent away from the solid, and viscous forces are reduced. 
If the plate is pulled, the interface is pulled toward the solid, making the film thinner and 
enhancing viscous effects. Hence in the latter case there is a positive feedback increasing 
viscous forcing and thus leading to instability. 

The strong energy dissipation near a moving contact line results from the fact that viscous 
forces become very large as the thickness of the liquid film goes to zero i^. As a result of 
the interplay between viscous and surface tension forces, the interface is highly curved, and 
the contact line speed U is properly measured by the capillary number Ca = Urj/'j, where 
rj is the viscosity of the fluid and 7 the surface tension between fluid and gas. Owing to 
this bending the interface angle measured at, say, 100/im away from the contact line differs 
[itI lisl . Iiol significantly from the microscopic angle directly at the contact line. 

The existence of a slip length A, which is of the order of a few molecular diameters under 

nnn 

normal circumstances j21l. 1221. 123|. implies a convenient separation of length scales between 
two different parts of the surface profile. On one hand, there is a contact line region whose 
typical scale is A, on the other hand there is an "outer" meniscus region, whose typical 
length scale is set by the capillary length £c = \J'^/{pg). Thus one expects the local "inner" 
behavior of the profile near the contact line to be of the form 

Mx)^3Ai^(f), (1) 

since A is the only available length scale. The dependence on the equilibrium contact angle 
6'p > was introduced for later convenience. 
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Firstly, one finds from (0) tliat tlie curvature of the interface is h"{x) = 6'^H"{$^)/{3X), 
which becomes large for A ^ as expected, while the slope is h'{x) = 6f.H'{C,). The inner 
solution ((H) has to be matched to an outer solution hout{x), whose curvature is of the order 
of which means we have to join the two solutions at some scale e with £c ^ e ^ A. This 
implies that the argument of hout can effectively be taken at x = 0, while e6'e/(3A) ^ 1. 
Thus the matching condition is 

h'U0)=0lH"{oc)/{3X), (2) 

which ensures that the inner and the outer solutions are compatible. We will see that ^ 
needs to be supplemented by what is essentially a condition for the slope. 

The matching condition Q is the key to understanding wetting behavior. Let us sum- 
marize the main results of this paper by analyzing the solution qualitatively for the two 
cases of an advancing contact angle (plate plunging into the fluid) and of a receding contact 
angle (plate being withdrawn). For very small A, it is clear that H"{oo) must be small for 
matching to be possible, so in the limit, H"{oo) = becomes the boundary condition for 
the inner problem. We will see that for an advancing contact angle this boundary condition 



yields the inner scaling function H{^), first found by Voinov 2M- The matching to Voinov's 
solution is slightly complicated by the presence of logarithmic terms in the slope jlfl| . 

On the other hand, Voinov's solution cannot be applied to the case of a receding contact 
line. Rather, all inner solutions maintain a finite curvature H"{oo) > 0. This means that 
at too small a value of A the matching condition (0) can no longer be obeyed and the inner 
and the outer solutions are incompatible. As a result, the contact line vanishes. In a typical 
experiment, A is constant and the speed is increased, but the effect is the same: since the 
curvature of the interface is caused by viscous forces H"{oo) increases with speed and 
matching becomes impossible. This explains the phenomenology described in Fig^ above. 
The impossibility of matching a receding contact line above a critical capillary number was 
already noticed in , based on numerical integration of the thin-film equations. Analytical 
solutions for the inner solution will permit us to give a much more complete description. 

In the next section we will introduce the hydrodynamic equations to be used for the cal- 
culation of stationary profiles. We will confine ourselves to the "lubrication approximation" , 
valid in the limit that the liquid film is thin. In the following section we consider the case of 
a solid plate being pushed into the liquid (advancing contact angle). By matching an inner 
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to an outer solution, we compute the profile as a function of speed. In the fourth section we 
introduce the matching procedure for the opposite case of a plate being withdrawn (receding 
contact angle). The failure of this matching procedure gives the critical capillary number 
at which the contact line can no longer exist, in agreement with our earlier result [l^. For 
speeds below the critical value we again find the interface profile. All our analytical results 
are tested by comparison with numerical solutions of the original equations. Finally, we 
summarize our results and indicate directions of future research. 



II. LUBRICATION DESCRIPTION 

As illustrated in Fig^ we are considering a plate being pushed into or being withdrawn 
from a liquid bath at an angle 9. This means we have to solve the steady Navier-Stokes 
equation with a free surface, and no-slip boundary conditions on the plate. Since the plate 
is moving with speed U, the contact line between fluid, solid, and gas is moving relative to 
the solid. As explained above, the Navier-Stokes equation does not allow such a solution, 
and a small-scale cutoff has to be introduced at the contact line jl?! 18^ . The dominant 



mechanism responsible for cutoff depends on the particular system under study |18l . |26| . As 
a representative example, we assume that the corner singularity is relieved by allowing the 
fluid to slip across the solid surface, since the matching to the slip region of size A is well 
understood in this case Ixo]. According to the Navier slip law |2], the fluid speed relative to 
the solid is proportional to the shear rate: 

du 

u{x, 0) - t/ = A— at y = 0. (3) 

Recently, we found ,27] that various modifications of Q have a minimal influence on the 
interface away from the contact line. Hence we do not believe that the particular cutoff 
mechanism used is of great importance. 

A much thornier issue is the slope h'{Q) of the fluid layer to be specified at the contact 
line. It is well appreciated that molecular processes are involved [3,13]' which are beyond 
a hydrodynamic description. This will lead to a an effective speed dependence of the contact 
angle as defined on a scale of nanometers. The importance of these microscopic effects 
relative to hydrodynamic ones is determined by the amount of energy dissipation involved 
in either process j30|. Thus for high viscosities the speed dependence of the microscopic 
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angle can most likely be ignored, in agreement with experimental data [31|. If there is no 
"intrinsic" speed dependence of the microscopic angle, it must coincide with its equilibrium 
value 6e at zero speed. We thus take /i'(0) = 6^. 

The mathematical problem simplifies significantly if we assume that the angle the inter- 
face makes with the solid is always small, and that viscosity is sufficiently large for inertia 
to be ignored. In this limit, one can find an approximate description of the hydrodynamic 



equations which eliminates the flow field 



ICt Il4| , and the so-called "lubrication equation" 



can be written entirely in terms of the free-surface profile h{x). Non-dimensionalizing all 
lengths with the capillary length £c = \/l / (pg) one finds [3] 

h'"-h' + 9, (4) 



/i2 + 3Xh 

where we consistently used the small-angle approximation tan(^) ^ 6. To distinguish more 
clearly between advancing and receding contact lines we always take Ca as a positive quan- 
tity, and rather change the sign in the equation. The - sign corresponds to the plate plunging 
into the liquid, the + sign to the opposite case of the plate being withdrawn. 

Viscous forces appear on the left of (jH) (proportional to the speed), and diverge quadrat- 
ically as h goes to zero at the contact line. As a result, viscous dissipation would diverge 
if it were not for the presence of slip, which weakens the singularity. Near the contact line 
the surface is highly curved, so the first term on the right of (jH), which comes from surface 
tension, balances the viscous term on the left. The other two terms stem from gravity and 
only come into play at greater distances from the contact line. The film thickness vanishes 
at the contact line, where the slope is h'{0) = Of. as discussed above. Far away from the plate 
the surface coincides with the liquid bath, so the third boundary condition is h'{oo) = 9. 
Rescaling the layer thickness h with the equilibrium contact angle Of.-, one finds that there 
remain three parameters in the problem, namely the combinations A/6'e, and d/9e. 

However, we find it more intuitive to keep the original parameters, and to state results in 
terms of the above combinations if convenient. 



III. PUSHING 



We begin by considering a plate being pushed into a viscous fluid, corresponding to the 
- sign in d^. In the spirit of the matching condition Q we approach this problem by 
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first considering the leading-order behavior near the contact hne, where h goes to zero. As 

ri 

discussed above, this equation is 

'^^"^ h'", (5) 



+ 3Xh 

dl in j^^l for a more general class of slip models. The characteristic 



which we studied in detail 
scale of the local solution is the slip length A, so it is convenient to introduce the scaled 
variables which leads to 



-H" 



(6) 



where 6 = 3Ca/6l is the rescaled capillary number. 

As argued above, the matching condition (0) leads to the requirement that the curvature 
H"{C,) vanishes for large ^ as the limit A — > is performed. This means that the boundary 
conditions for the solution of © are 



H{0) = 0, H'{0) = 1, H"{oo) 



0. 



(7) 



The only parameter now appearing in the problem is the rescaled capillary number S, and 
equations (lEI)-© uniquely specify the profile close to the contact line. 

This inner solution can be found bv exp an ding in a power series in the capillary number 
in a manner described in many papers |l0l.l27ll32j. If one writes the solution in terms of h'^, 
its behavior for large x/X can be written as 



Kn{.x)-ei = 9Caln{x/L) 



(8) 



io any order in the capillary number [27]. This is the form originally proposed by Voinov 



2J| , using more qualitative arguments. The length L appearing inside the logarithm can be 



computed as a power series in the capillary number: 

^2 1 



TT 



2^3 



Ca + 0{Ca' 



(9) 



but for simplicity we only take the leading order term into account here, as corrections 
introduced by the next order are usually quite small [27 1. 

The crucial point of representing the inner solution in the form (jH)) is that the only 
parameter multiplying ln(x) is the capillary number, which is defined in terms of the "outer" 
problem. Namely, the outer problem (jH)) needs to be matched to is 

-3Ca 



"-out 



h 



out 



Kut + ^5 



(10) 
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which does not contain any contact hne parameters hke A or 9^. Owing to the strong 
singularity for — > [3, ()10p does not have a solution with finite slope at the contact 
line hout{0) = 0, making it impossible to impose a slope. Instead, we are seeking a solution 
of that has the form (jH} for x — 0, in other words 

h'^^tix) = 9Ca ln{x) + F, x ^ 0, (11) 

where F is a constant to be computed. As demonstrated in [l^, this can be achieved by 
expanding the outer solution in a power series in Ca: 

houtix) = ho{x) + Cahi{x) + 0{Ca^). (12) 

The equation for ho, representing a balance of surface tension and gravity, is 

h'^ -h'o + e = 0. (13) 

Equation (fTT^ has the following family of solutions, which are finite at infinity and which 
vanish at the contact line: 

ho{x) = ex + {e-eap){e-^-i). (w) 

The slope h'^uti^) = ^ap of at the contact line is called the "apparent" contact angle and 
is a free parameter still to be determined by the matching procedure. Its name is motivated 
by the fact that a macroscopic measurement of the profile on the scale ic will yield a profile 
close to ()14|) . If extrapolated to the contact line, the angle will appear to be 6ap, rather than 
the true microscopic value 6e- Obviously, p4|) cannot be matched directly to the contact 
line solution (jH)), which contains a logarithm. This logarithmic dependence will come out of 
the next-order solution hi, whose equation reads 

h'^-h\ = f{x), fix) = ~3/hoixf. (15) 

The solution is straightforward jSSjJ: 

hi{x) = -^J e-'f{t)dt + J {cosh{t - x) - 1) f {t)dt + K e-"" + K2, (16) 

where once more any particular solution that is growing at infinity was suppressed, and K, K2 
are constants of integration. To find K, K2, we note that h must vanish at the contact line, 
giving the first condition hi{0) = 0. To compute /ii(0), we observe that the second integrand 



9 



of ()16p can be expanded like cosh(t — x) — 1 = cosh(t) — 1 — x sinh(t) +0(x^). The term linear 
in X does not contribute for small x, since sinh(t)/(t) behaves like 1/t for small arguments, 
and hence 

/•X 

limx / smh{t)f{t)dt = 0. 

Thus one finds 

-| roc /»1 

= /ii(0) = — / e^'f{t)dt+ {l-cosh{t))f{t)dt + K + K2, (17) 

where the singularity of f{t) oc cancels out to make the second integral convergent. 
Next we find a condition at infinity by noting that 

h{x) = 9x + A + 0{l/x), (18) 

where A is vertical distance of the contact line from the undisturbed surface. Since we want 
ho to represent the far-field behavior of the profile, we put hi{oo) = 0, and thus A = 9ap — 0, 
by comparison with ()14|1 . This means we have 

= /ii(oo) = Jim e*/(t)rft|-^ f{t)dt + K2. 

The first of the two terms on the right is zero, as one confirms by splitting it into two parts: 

-.-A rB ^-A rA 



Jim e'f{t)dt+^— e'f{t)dt^ 



where i? is a large positive constant. The first of the two parts is evidently zero, for the 
second one notes that the argument can be approximated as e^f{t) ~ —3e*/{6'^t'^) for large 
t > B, whose absolute value has the upper bound 3e^ / {O"^ A^) . Thus in the limit the second 
part vanishes as well. Using ()17p this means the constant K in ()16p can be computed as 

/oo -t rl 
(— - l)f{t)dt + (cosh(t) - l)f{t)dt. (19) 

We are now in a position to determine the constant F in (fTT|) . Comparing (fTT|) to (fT^ . 
we know that for small x the first order contribution h'^ must have a logarithmic singularity 
of the form 

K{x) ^ 3/el^Hx) + c{e/e,,)/el, (20) 

Namely, for small x 

9Caln(x) h'^^^{x) ^ {h'^^x) + Cah[{x)f ^ 6% + 9Caln(a;) + 3CaC(^/0ap) + 0{Ca^), 
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FIG. 2: The depression A of the meniscus (divided by the equihbrium contact angle 9e), as 
function of the reduced capillary number 3Ca/6^. The slip parameter is 2>\/9e = 10~^, 6/9(, = 1. 
The full line is the result of a numerical solution of (01), the dashed line is our theoretical result 

n 

((221); using A = Oap — 9. The dotted line is the theoretical result of 

so that F = dlp + 3CaC{9/9ap). Analysis of dEl) for a; gives, using ^ and f{t) ^ 
-^/(dapt)^ for small t, 



0(9/9. 



ap J 



(1 - e'')fit)9ldt + [(1 - e-')f{t)9l + 3/t] dt + 3. (21) 



This completely determines the the outer solution (fT^ . Comparing it to the inner solution 
(jH)) we finally obtain 

9%^ + ?>CaC{9/9ap) = 91 - 9Caln(L), (22) 



which is an equation to be solved for the apparent contact angle 9ap- 

To test the result of our matching procedure ()22j) . we compare the depression A = 
9ap — 9 of the meniscus with the result of a numerical solution of the original equation (^. 
Remarkably, the prediction, which contains no adjustable parameters, remains extremely 
good up to a reduced capillary number of 3Ca/9l = 1. To further appreciate the quality of 
the agreement, we plotted the result of an earlier theory ^| as the dotted line. In this earlier 



11 




4e+05 



xe^/3?i 



FIG. 3: A comparison between a profile obtained numerically by integrating Q) at 5 = ?>Ca/6l = 
0.3 and the outer solution H12|) . The other parameters are those of Fig|2 For 5 = 0.3, (|22() gives 

Oap/ Ge = 2.32. 

theory, a perturbation expansion of (jH) in Ca is matched directly to a linearized version of 
(jH)). However, the expansion is performed around the static profile corresponding to zero 
speed, whereas our expansion is around ho^x), which already incorporates a speed-dependent 
deformation of the surface. For the special case 6 = 6e, the result of ll| is 



A = -9 J ML) - 7] 



(23) 



where 7 is Euler's constant. A comparison of the profile h'{x) obtained numerically with the 
outer solution (fT^ shows equally good agreement, as shown in FiglSl It is only for x/A < 1 
that the outer solution begins to fail, since it has a logarithmic singularity at the origin. 
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IV. PULLING 

We now turn to the opposite case of a plate being withdrawn from a bath, for which the 
left-hand side of dH) is now positive. The dominant balance close to the contact line is again 
between the term on the left of ^ and the first term on the right. With the similarity 
transformation ((T)), this is converted into 

which differs from only by a sign. However, the behavior of solutions of (j24j) as ^ ^ oo 
is completely different. This is best appreciated by considering the form of (j24p valid away 
from the contact line where H ^ 1: 

= y'", (25) 

where we have put H{^) = 5^^^y{^). 

Jlemarkably, this equation has an exact solution, whose properties have been summarized 
in In parametric form, a solution with y{0) = reads 



y 



l3{aAi{s)+l3Bi{s)) 
1 



s e [sicx)[, (26) 

where Ai and Bi are Airy functions 37]. The limit ^ — > corresponds to s — > oo, the 
opposite limit ^ ^ oo to s si, where si is a root of the denominator of (|^: 

aAi{si) + PBi{si) = 0. (27) 

Since the solution extends to s = cxd, si has to be the largest root of (j77|) . 

From the behavior of y (^) for large ^ can be obtained 3^ (note that there is a 
misprint in equation (12) of jSJ]): 

y'{0 = '^y^ + by + O{r'), (28) 

where 



" \t:Ai{s{)) ' " Ai{s^) ■ 
The constant (3 can be determined by matching (pUj) . which is valid only for ,^ > 1, to a 

solution of TIM , which includes the effect of the cutoff and is thus valid down to the position 

n 

^ = of the contact line. The limit of for small values of ^ gives j2J| 

H'\x) = 5y'\i) ^ 351n(7r/(22/3/52^), (29) 
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which remains vahd for ^ 13 Thus, ()29|) is a vahd solution of the full equation ()24p for 

1 < ^ < 

Following |l0l|27|, we compare (|29|) to the expansion of the full equation (|24j) in 6. Using 



the boundary conditions H{Q) = and H'{0) = 1, one finds 

H'io = 1 + '5 laHO - in(e + 1)) - in(e + 1) + ce] + o{s'). (30) 

Since this solution has to match the logarithmic behavior (j^^ . we put the constant of 
integration C to zero, and (jHUj) becomes for ,^ ^ 1 

i/'3(0 = l-351n(0 + O(52). (31) 

Thus, comparing (|!TT|l to we find 

= 7rexp (-1/(35)) 72^/3 + 0(5). (32) 



For S <^ 1, P is indeed exponentially small, and ()26|) has a logarithmic dependence over a 
large range of .^-values: 1 ^ ^ 6xp(l/(35)). This makes it possible to match to (jHTj) 
in the limit ^ ^ 1, although the ultimate behavior of ()26|) for large ^ is given by ()28|). 

We are now in the position to match the inner solution of ()24j) . which has the form (P), 
to an appropriate outer solution, following the prescription Q given in the introduction. To 
achieve this matching it is enough to consider hout{x) = ho{x) as given by (fT^ . since 
does not contain any logarithmic term. For x/X ^ 1, the solution of is well described 
by (j2ni), so we have 

houtix) = 9x + {9 - 9ap) 1) (33a) 
h^x) = 3X6'/'y{x9J{3X)), (33b) 

where y{^) is given by (j^ . ljH^ . The matching procedure must supply us with the parameter 
si in p6|l . which is yet to be determined. 

The matching works in the limit A 0, for which we require that the expansions 

KUx) = 9,, + {9-9,,)x + 0{x^) 
h[^{x) = 9J'/^ [Kyx9J{3X) + by] + OiX/x) 

agree, or 



(34) 



9eS^%, (35a) 
9-9ap = 9y/\y/{3X). (35b) 
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FIG. 4: The apparent contact angle 9ap according to pZf) . The other parameters are \/9e = 10~^ 
and 0/9(. = 2. The dashed curve is the approximation ()44() . vaUd for smah capillary number. The 
apparent contact angle goes to zero at the critical capillary number Cocr/^f = 0.0103. 



Eliminating 9ap between the two equations we finally have 

6 22/3Ai'(si) ^eexp(-0,V(9Ca)) 



(36) 



(3Ca)V3 ^^(s^) 3-2i/37rA?2(si)A ' 

which should be read as an equation for Si. Once Si is known, all remaining parameters in 
hin and hout can be found. For example, from ()35a|) one finds the apparent contact angle to 
be 



e, 



ap 



-2^/\'iCaYl^Ai\s^) 



(37) 



an example of which is plotted in Fig. |31 At some finite critical capillary number Cacr, 
6ap goes to zero. Physical solutions cannot exist for capillary numbers beyond that, since the 
outer solution only makes physical sense for 9ap > 0. Following (l3|, we now derive a much 
simpler equation for the critical capillary number Cacr at which the contact line disappears, 
and show that it is equivalent to 6ap = 0. 

Since 9ap > 0, the maximum value of the left hand side of ()35bj) . which is the curvature 
of the outer solution, is 6. The right hand side, corresponding to the curvature of the inner 
solution, is on the other hand bounded from below. From ()28|) it is seen that the minimum 
curvature corresponds to the condition that Ai[si) must be maximal among solutions of 
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(j2Zj). By choosing a = = —PBi{smax) / Ai{smax) we can in fact ensure that Ai assumes 
its global maximum 0.53566 . . . , which occurs for s = Smax = — 1-0188 .... Thus we have 
singled out a unique solution of which minimizes the curvature 

^ exp[-gf/(9Ca)] 
''y 2yMMsmax)y' ^ ^ 

the value of which increases with capillary number as expected. 

Now by equating the maximum value of the left hand side of (j35bj) with the minimum 
value of the right of ()35bj) we obtain an equation for a capillary number above which no 
solution can exist: 9 = 916^/^kI''/{3X), or exphcitly 
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In 



C (Xcr 



-1 



(39) 



But at the capillary number given by ()39|1 . the first matching condition ()35aj) is also satisfied 
identically, since 6ap = by = 0. This is because Ai is extremal, so Ai' = 0. Thus Cacr as 
given by (jH^ gives exactly the critical capillary number corresponding to 6ap = in Figl^ 



This confirms a classical conjecture by Derjaguin and Levi |35|, later reiterated by others 
j^^, that the transition to a film is characterized by the apparent contact angle going to 
zero. This criterion was confirmed experimentally in jj^, using fibers being pulled out of a 
viscous liquid. 

The concept underlying equation ()39|) for the critical capillary number is illustrated again 
in FigI3 Below the critical capillary number, there is some overlap in the curvature of the 
inner and the outer solutions (left diagram in Fig|S)). The matching equation (jHUj) determines 
which solution is selected from this overlap region. Above the critical capillary number, 
there can be no solution (right diagram in FigEl): there is no outer solution which could be 
matched to an inner solution, because this would result in a sudden jump in the curvature. 

In Fig. El we show the result of a numerical integration of (jU at the critical capillary, 
and compare it to the inner and outer solutions (jHSJ- The critical capillary number was 
found numerically by raising Ca until no more solutions of Q could be found. For large 
^, the slope of the outer solution ()33a|l agrees with the numerical solution, but extrapolates 



to h' = at the contact line, as reqired by 6ap = 0. Coming from the interior, the inner 



solution (j33b|l agrees with the full solution up to the turning point. 

To obtain a profile that is more uniformly valid, one can use a composite approximation 
j^. The idea is to add the inner and outer solutions that have been matched, and to subtract 
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FIG. 5: A schematic illustrating the transition to a LLD film for a receding contact line. Below 
Cccr, the curvatures of all possible inner and outer solutions have an overlap region, so they can 
be matched. Which solution from the overlap region is selected (illustrated by the dotted line), is 
determined by (|36() . Above Cccr no matching is possible, since the the maximum curvature of all 
possible outer solutions lies below the minimum curvature of all the inner solutions. 

the behavior (|34j) in the region in which the two solutions overlap: 



It is evident that hcomp agrees with the full solution for large as well as for small x, hence 
it will be the best global approximation at this order of the matching. This excludes the 
region ^ < 1 where ()33b|) has a logarithmic singularity. Figure [7| illustrates the remarkable 
agreement of the composite solution with the full numerical result at the critical capillary 
number, where the approximation is expected to be worst. 

Having studied the critical capillary number, let us return once more to the matching 
condition ()3(ij) which determines the solutions for Ca < Cacr- Namely, it is instructive to 
obtain explicit solutions of in the limit of small Ca, in which case its solution Si is 



hcompix) = hin{x) + hout{x) - [9ap + (9 - 9ap)x] . 



(40) 
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FIG. 6: A comparison of the full solution at the critical capillary number with the inner and outer 
solutions. We plot the slope of the interface, so h'{Q)/9e = 1 for the full solution, and /i'(0) = 
for the outer solution, consistent with the condition by Derjaguin and Levi. The other parameters 
are X/Oe = lO'^ and e/Oe = 2, so from ^ we obtain Cacr/Ol = 0.01032. 



large, and the asymptotics of the Airy function j37| gives 

Ai{s)^e-^'"'"l^s-^'^/{2'R^/^). (41) 
Thus, (jSni) becomes for small Ca: 

and (jHTf) is 

e^p ~ (3Ca)i/322/3^i/2_ (43) 

In the simplest approximation the exponent in (jl^ must vanish, s = {6l/{12Ca))'^^^ , and 
inserting this into (j43|) gives O^p ~ Og. Not surprisingly, at very small capillary numbers the 
interface is hardly deformed at all. 

To go beyond this approxiamtion, we put s^^^ = (^g/(12Ca) + As)^^^ , insert into ()42p . 
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FIG. 7: The interface slope at the critical capillary number for 0/6e = 2 and \/6e = -10"^. The full 
line is the result of the numerical integration at Cacr/Ol^ = 0.0103, while the theoretical prediction 
gives Cacr/6^ = 0.01032. The dahes line is the composite solution H40() . 

and analyze the result for small Ca, which gives 

As - In 

4 

Putting this into ()4H|1 leads to the next approximation 

which is similar to equation ()22|) for the apparent contact angle in the advancing case. 
Indeed, for small capillary number ()22|) becomes 

which has the same form as (j44j) . However, even for small capillary number the receding 
and the advancing case needs to be treated differently, as the constant inside the logarithm 
is different. For 9 = 9^ the argument of the logarithm in pij) becomes singular, which 




3A 
2^' 



29i 



3A 



(44) 



egeexp(C(g/ge)/3) 
3A 



(45) 
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means one has to go to an even higher order in the approximation, but we are not pursuing 
this special case here. The approximation ()44|) is plotted as the dashed line in FigEJ If 
extrapolated naively to 9ap = 0, it gives a surprisingly good estimate of Cacr- 



V. DISCUSSION 



Let us begin by considering rather straightforward extensions and generalizations of the 
present theory. Firstly, our arguments are not limited to a specific contact line model, 
since they are based entirely on hydrodynamic arguments away from the contact line. For 
example, if van-der-Waals forces are dominant near the contact line this only changes 
the parameter L appearing in (jH)) (advancing contact angle) or f3 in ()29p (receding contact 



line). In the latter case, the slip length A in (jH^ has to be replaced by A/{6TT'y)/{Q9e), 
where A is the Hamaker constant. Our research also suggests that the calculation of the 
critical capillary number is captured fully by lubrication theory, even if 6e is not small. In 
that case the flow directly at the contact line would have to be described without resorting 
to lubrication theory Q^. but the relevant region where matching occurs . charaoter.zed 
by slopes h'{x) which are small. 

Secondly, one can generalize to a different geometry. To this end one has to replace ()14|) 
by the appropriate static solution for the problem at hand. This is done almost trivially 
for the case of a vertical plate or a fiber , in which case the lubrication description (j3)) is 
no longer valid far from the contact line, since 6 is not small. This however does not pose 
a problem since this part of the profile is determined by surface tension and gravity alone. 
In the same spirit, the present model can be exteded to a flow inside a capillary tube, with 
only very minor changes to the value of the critical capillary number Q|. 

Thirdly, one can consider dvnamical effects, of particular interest for the unstable case 
of a receding contact line. In [lj| it was found that when the stationary profile vanishes, it 
is not followed directly by the LLD film, which is of macroscopic thickness. Rather, there 
is a narrow range of speeds where the contact line is pulled up the plate, but at a speed 
that is smaller than U, i.e. the contact line is partially slipping. The thickness of the film 
that is left behind is in the order of A, i.e. microscopic. Only when the speed is raised still 
further does the LLD film appear. So far these results are only numerical; a full analytical 
theory would be desirable. Note also that these transitions are strongly hysteretical. Once 
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the LLD film has appeared, it can be sustained to much lower capillary numbers than Ca^r- 
It is usually assumed that the LLD film vanishes when its thickness has reached the range 
of intermolecular forces jj^], but we are not aware of any theoretical investigation of this 
problem. 

Next we come to the experimental evidence. The most extensive experiments on the 
critical capillary number for a receding contact line were done with a capillary tube Q|, 
by pushing out a viscous liquid. Using a variety of different materials it was found that 
instability occurs at a given value of the reduced capillary number Ca/6'g, in agreement 
with the present theory. The actual value of the critical capillary number, however, is 
about a factor of two too low, if the theoretical estimates are based on A ~ nm. Several 
pos.b.e e.p>a.at,o,. suggest tKe,nse.ves. Fi.t, the .natetials used ,n Q have eous.de.ab.e 
contact angle hysteres is, p ointing to surface roughness. This will tend to reduce the critical 
capillary number js^ [4^. In addition, any speed dependence of the microscopic contact 
angle, neglected in the present description, will effectively lower Of. and thus lead to a smaller 
critical capillary number. 

Our theory for the vanishing of the receding contact line has some similarities with an 



earlier theory 



4l|, in that the critical capillary number is proportional to 9^ (cf. (jSHl))- 



However, it differs in predicting a vanishing apparent contact angle at the transition, while 
it is Bap /Of. = l/VS in |4]J. The approach of 4l| is also different in that it considers the 



local problem in isolation, hence the dependence of Cacr on parameters of the outer problem 
like 6 cannot be captured. In fact, we believe that the mechanism for instabilitv proposed 
, which is based on an approximate solution of (j23), contains a flaw. In 27] we use 
the case of the advancing contact angle (0) to show that the method of solution proposed 



m 



4l| cannot correctly predict the nonlinear dependence of the angle on speed. But it is 



4l|. 



precisely this nonlinear dependence which lies at the heart of the stability analysis of 

The most important extension of the present theory of contact line instability however 
is its application to higher dimensions, in which the contact line no longer remains straight. 
If the plate withdrawn from the liquid bath is sufficiently wide, the contact line inclines 
relative to its direction of motion. Two sections of the contact line that have inclined in 
opposite senses meet at a sharp corner, so that the whole contact line is serrated in an 
irregular fashion js^. 

A more controlled recent experiment is that of a viscous drop running down an inclined 
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plane [3|. At a critical speed, the initially rounded tail of the drop forms a sharp corner. 
A recent theory explains the drop profile near the corner of the drop, but not the 
critical speed at which the the corner first occurs, nor its opening angle. To give a complete 
description of the transition, the microscopic neighborhood of the contact line has to be 
included, as was done in the present theory for a straight contact line. One important 
difference between the case of a sliding drop and that of the present paper is that recent 
experimental evidence suggests that the transition occurs at a finite value of the apparent 
contact angle Studies to understand the three-dimensional instability of a sliding drop 
are currently under way. 
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